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Summary 

Low-frequency  attenuation  in  the  Arctic  Ocean  is  more  than  an  order  of 
magnitude  greater  than  that  in  the  open  ocean.  Scattering  at  the  ice/water 
interface  is  the  most  likely  mechanism.  Earlier  estimates,  based  on  Marsh 
theory  for  sea-surface  scattering,  showed  reasonable  consistency  with 
experiment.  However,  analysis  of  more  recent  experiments,  with  roughness 
parameters  derived  from  concurrent  measurement  of  the  underice  profiles, 
yields  values  that  are  low  by  more  than  a  factor  of  two. 

The  approximations  of  Marsh  theory  have  been  investigated.  More  exact 
calculation  of  the  scattering  integral  by  numerical  methods  has  been  found 
to  increase  the  discrepency  at  low  frequencies. 

Failure  of  the  pressure-release  approximation  to  account  for  attenuation 
suggests  that  impedance  effects  are  involved.  For  impedance  calculations, 
the  ice  layer  is  modeled  as  a  uniform  elastic  solid,  assuming  that  acoustic 
contact  at  the  ice/water  interface  is  perfect.  Trapped  air  could  reduce  the 
effective  impedance,  making  the  pressure-release  assumption  valid.  The 
transfer  function  of  acoustic  pressure  to  vertical  velocity  is  estimated  by 
comparing  hydrophone  and  seismometer  signals  and  no  evidence  of  such  an 
impedance  discontinuity  is  found.  However,  horizontal  surface  velocities 
are  much  greater  than  predicted  by  thin-plate  theory,  which  supports  a 
"rocking"  motion  produced  by  slopes. 

Perturbation  analysis  shows  scatter  by  both  the  vertical  displacements 
and  the  slopes  of  the  ice/water  interface.  The  displacements  produce  a 
vertical  velocity-dipole  component  of  scatter  in  addition  to  a  vertical 
pressure-dipole  component  analogous  to  the  free-surface  case.  The  latter 
is  the  dominant  mechanism  at  higher  frequencies  while  the  former  is 
negligible.  The  horizontal-dipole  component  produced  by  slopes  may  be  a 
significant  mechanism  at  low  frequencies;  however,  the  predicted  values 
of  scattering  loss  still  fall  far  short  of  agreement  with  data. 
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Introduction 


Sound  Speed 
(km/s) 


Figure  1:  Arctic  Ocean  sound-speed  profile  and  ray  diagram. 


The  Arctic  Ocean  is  a  unique  acoustic  environment,  principally  because 
the  sound  speed  increases  monotonically  with  depth  and  the  surface  is 
usually  covered  by  ice.  Upward  refraction  coupled  with  the  effects  of  the 
rough  ice/water  interface  cause  severe  reverberation  and  attenuation, 
both  of  which  can  limit  SONAR  performance. 

Figure  1  shows  a  typical  sound-speed  profile  and  ray  diagram.  In  deeper 
waters,  the  maximum  RSR  cycle  distance  is  of  the  order  of  50  km.  At  short 
range,  reverberation  and  absorption  in  the  medium  are  the  limiting  factors. 
At  long  ranges,  most  of  the  energy  is  concentrated  at  low  grazing  angles. 
Multiple  reflections  at  the  ice/water  interface  are  then  involved  and  the 
resulting  loss  can  exceed  absorption  by  more  than  an  order  of  magnitude. 

The  suspected  mechanism  is  scattering  at  the  rough  ice/water  interface. 
Early  modeling  efforts  were  based  on  the  sea-surface  theory  of  Marsh  111. 
Free-surface  theory  appeared  to  give  plausible  results;  i.e.  the  predictions 
were  consistent  with  experiment  for  the  standard  deviation  2.3m,  a  value 
in  reasonable  agreement  with  underice-profile  data. 

Marsh  theory  is  valid  only  over  a  limited  range  of  angle  and  frequency. 
Recent  advances  in  scattering  theory  now  permit  numerical  calculation 
over  a  wider  domain.  In  addition  to  improved  theory,  data  are  available 
from  underice  profiles  taken  concurrently  with  propagation  measurements. 
With  more  reliable  estimates  of  standard  deviation  and  correlation  length, 
the  calculated  scattering  losses  for  the  free-surface  approximation  are 
consistently  low  by  more  than  a  factor  of  two. 


In  earlier  work,  the  ice/weter  interface  was  taken  as  compliant  mainly 
because  theory  for  finite-impedance  boundaries  was  not  available.  Theory 
for  a  rough  fluid/solid  interface  still  does  not  exist.  To  approximate  the 
effects  of  acoustic  impedance,  we  consider  a  uniform  flat  elastic  layer 
and  treat  roughness  by  the  perturbation  method.  Impedance  is  calculated 
by  the  theory  of  Brekhovskikh  [2], 

In  the  compliant  case,  displacements  reradiate  as  vertical-dipoles.  At 
low  frequencies  and  grazing  angles,  ice  impedance  tends  to  mass-loading 
and  the  elastic  effects  vanish.  Ice  floats  on  water  and  is  isostatic  on  the 
average,  the  total  mass  being  equal  to  that  of  the  displaced  water.  In  this 
limit,  the  layer  would  become  quite  transparent.  However,  if  isostacy  does 
not  occur  on  the  scales  involved,  some  vertical-dipole  scattering  can  be 
expected  since  the  ice/water  density  ratio  is  roughly  0.9. 

Although  the  roughness  of  the  air/ice  interface  can  also  have  an  effect, 
the  magnitudes  are  small  because  standard  deviations  are  small  compared 
to  the  ice/water  interface. 

Losses  due  to  dissipation  within  the  layer  also  appear  to  be  small  except 
at  high  frequencies  and  grazing  angles  where  resonance  occurs  between 
the  acoustic  wave  in  water  and  the  shear  or  compressional  waves  in  ice. 
For  the  low  frequencies  and  grazing  angles  involved  in  Arctic  propagation, 
dissipation  should  be  negligible. 

Impedance  effects  could  become  negligible  if  acoustic  coupling  at  the 
ice/water  interface  were  reduced  by  trapped  air.  Comparison  of  vertical- 
seismometer  and  hydrophone  signals  shows  no  indication  of  this  type  of 
impedance  discontinuity.  However,  the  horizontal-seismometer  signals  are 
much  greater  than  predicted  by  flat-plate  theory,  indicating  that  slopes 
are  also  involved  in  the  ice  motion  in  response  to  acoustic  waves. 

The  perturbation  method  does  not  consider  local  grazing-angle  effects. 
Finite  slopes  increase  the  "effective”  grazing-angle,  which  also  increases 
acoustic  coupling  at  rough  solid  boundaries.  Schmidt  and  Kuperman  (3)  have 
considered  this  problem;  however,  their  analysis  is  based  on  the  Kirchhoff 
(tangent-plane)  method,  which  is  valid  only  at  higher  frequencies. 

At  finite-impedance  boundaries,  the  slopes  produce  a  horizontal-dipole 
component  of  radiation  associated  with  "rocking"  motion.  Analysis  shows 
that  this  is  a  significant  scattering  mechanism  at  the  lower  frequencies. 
However,  attenuation  calculations  still  fall  far  short  of  agreement  with 
experimental  data  in  this  range. 

The  purpose  of  this  report  is  to  review  the  experimental  and  theoretical 
evidence,  in  the  hope  that  resolution  of  the  remaining  questions  will  lead 
to  an  accurate  predictive  model. 


-  2  - 


Attenuation  data 


Figure  2:  Arctic  attenuation  measurements. 

figure  2  compares  attenuation  data  with  predicted  absorption  in  Arctic 
sea-water.  The  solid  circles  are  1959  NUSC  dat8  [4]  and  open  circles  are 
Tristen/Frarn-82  data  151.  The  connecting  curves  are  arbitrary. 

The  striking  features  of  the  dat8  are  the  magnitude  and  resemblence  to 
relaxational  phenomena.  Below  1  kHz,  the  values  are  more  than  10  times 
greater  than  absorption.  Saturation  at  high  frequencies  is  believed  to  be 
due,  in  part,  to  experimental  conditions.  At  high  frequencies,  attenuation 
limits  the  useable  range  to  less  than  50  km  and  absorption  dominates  in 
this  regime.  At  low  frequencies,  long  ranges  are  required  in  order  to  make 
accurate  measurements  and  ice  loss  then  becomes  dominant. 

The  likely  mechanism  Is  scattering  at  the  rough  ice/water  interface.  The 
theory  for  free-surface  scattering  shows  that  the  reflection  loss  depends 
linearly  on  grazing  angle  for  small  angles.  Since  skip  distance  has  similar 
dependence,  loss  vs  range  depends  mainly  on  gradient.  Ray  angles  greater 
than  about  8°  encounter  a  smaller  gradient;  however,  divergence  effects 
are  greater  and  contributions  are  small  at  long  ranges.  Cutoff  by  bottom 
topography  may  be  an  important  factor. 

We  next  examine  the  statistical  model  of  underice  roughness  to  be  used 
in  calculating  reflection  loss. 
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9  9  2  3/2 

1- Dimensional  Spectrum  S,  =  li  qJ/(qJ+q  ) 

Correlation  Function  ^  =  h2(Qor^i 

2- Dimensionol  Spectrum  S2=  2  h2fto/(ftQ+ft2)2 

Standard  Deviation  RMS  h=2m 
Correlation  Length  L=2/Q0=44m 


Figure  3:  Ice  roughness  spectrum  model. 
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Figure  3  shows  o  sketch  of  a  typical  ice  profile.  The  graph  is  a  spectrum 
of  the  ice/water  interface  roughness  (Air/ice  interface  roughness  is  small 
and  will  be  neglected).  The  analytic  spectrum  (smooth  curve)  is  derived  by 
curve-fit  of  experimental  spectra  and  is  given  by  St  =hr2 fJo2/(  Q2  +  Q2^2 

where  Q  is  the  angular  spatial  frequency  or  wavenumber.  The  correlation 
function,  derived  by  cosine  transform  of  5V  is  'F(r)=h2(f2or)  Kt(n  r),  where 
K(  is  the  modified  Bessel  function,  L=2/Qo  is  the  correlation  length  and  h 

is  the  standard  deviation  about  mean  draft.  Azimuthal  distribution  is  taken 
to  be  isotropic  and  the  2.-dimensional  spectrum,  derived  by  the  Jo  Bessel 

transf orm,  becomes  S~=2  h2  Q- 2/(  Q  2  +  ft2)2  where  [°°Q  dQ  S0(Q)=h2 

Ice  profiles  of  draft  were  analyzed  in  approximately  1  km  segments  and 
averaged  over  10  contiguous  segments.  The  power  spectra  and  correlation 
functions  were  calculated  numerically.  Correlation  lengths  were  measured 
by  fitting  correlation  curves  with  the  analytic  model. 


Mean  Draft  (m) 


Figure  4:  Standard  deviation  vs  mean  draft 

The  standard  deviation  vs  mean  draft  data  in  Figure  4  show  quite  high 
correlation.  The  RMS  value  for  the  set  is  h=2m.  Figure  5  shows  negligible 
correlation  between  L  and  mean  draft.  The  average  of  the  set  is  L=  44m. 
Figure  6  shows  no  significant  correlation  between  L  and  h  also. 


Mean  Draft  (m) 


Figure  5:  Correlation  length  vs  mean  draft 


Figure  6:  Correlation  length  vs  standard  deviation 
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Figure  7:  PDF  of  draft-model  and  data 

Analysis  of  a  "shot-noise"  model  of  the  underice  roughness  by  Nuttall  [6] 
indicates  that  the  probability  density  function  (PDF)  of  dr8ft'  should  be 
similar  to  a  Rayleigh  distribution. 

In  Figure  7,  the  expected  PDF  (smooth  curve)  is  compared  with  typical 
data.  Both  have  been  normalized  to  unit  area  and  have  the  same  standard 
deviation.  Note  that  the  mean  draft  of  the  model  is  clearly  greater  than 
that  of  ice  and  the  PDF  is  less  highly  peaked.  The  purpose  of  the  composite 
model  is  to  reconcile  these  discrepencies. 

The  high  degree  of  correlation  between  standard  deviation  and  mean 
draft  in  Figure  4  indicates  that  the  underice  roughness  is  not  statistically 
homogeneous.  There  appear  to  be  significant  variations  in  the  roughness 
parameters  that  occur  on  a  scale  much  larger  than  the  correlation  length 
and  this  suggests  that  modeling  roughness  as  a  single-population  may  not 
be  appropriate. 

The  ice  pack  shows  large  areas  of  comparative  smoothness  interspersed 
by  areas  of  high  roughness.  We  will  consider  a  composite  model  composed 
of  two  populations  having  different  degrees  of  roughness  and  occupying 
different  contiguous  areas.  Note  that  this  differs  from  the  conventional 
sea-surface  composite  model  in  which  two  are  superposed. 
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Figure  6:  PDF  of  draft-composite  model  and  ice  data 

In  Figure  8,  the  data  are  fitted  using  two  distributions  having  different 
standard  deviations,  i.e.  1.1m  and  3.3m,  shown  by  the  dashed  curves.  The 
weighting  factors  are  0.7  and  0.3,  respectively.  A  bias  of  1.5m  has  been 
added  to  improve  the  fit  and  the  respective  mean  drafts  are  approximately 
4m  and  8m.  The  sum  of  the  two  components,  shown  by  the  solid  curve,  is 
now  in  fairly  good  agreement  with  the  data. 

Scattering  loss  is  proportional  to  the  variance  h2.  When  the  correlation 
lengths  of  the  two  components  are  equal,  the  area-weighted  variances  are 
additive.  The  variance  of  mean  draft  has  a  much  greater  correlation  length, 
which  will  reduce  its  contribution  in  Marsh  theory.  This  component  will 
appear  as  an  increase  in  spectrum  level  at  very  low  wavenumbers  and  a 
slow  decay  in  the  correlation  function  at  long  lags.  Both  effects  have  been 
noted  in  the  experimental  data. 

While  the  composite  model  appears  somewhat  more  realistic  than  the 
homogeneous  model,  the  effects  on  scattering  appear  to  be  negligible,  8t 
least  in  the  free-surface  case.  Concentration  of  roughness  in  one  area  at 
the  expense  of  another  does  not  increase  low-frequency  loss  significantly 
in  this  approximation.  It  will  be  seen  later,  however,  that  the  effects  do 
become  important  when  impedance  is  taken  into  account. 


-  8  - 


Scottehng-loss  theory 


The  Marsh  formula  for  the  coherent  energy  reflection  coefficient  can  be 
written  as  Q=1-SL,  where  SL-3.3  sin8  (kh)2/(kL)1/2,  8  is  grazing  angle  and 
k  is  acoustic  wavenumber.  To  allow  some  extension  beyond  the  limits  of 
perturbation  theory,  the  loss  is  approximated  as  4.34  SL  (dB/bounce). 

The  skip-distance  in  a  surface-channel  is  given  by  2c  tan8/g  where  c  is 
the  sound-speed  and  g  is  the  gradient.  Loss  vs  range  is  then  approximately 
independent  of  ray-8ngle  for  small  angles.  For  g=.06/sec,  h=2m  and  L=44m, 
the  attenuation  becomes  As  1.5  f3/2  dB/km  where  f  is  the  frequency  in  kHz. 
This  result  is  in  good  agreement  with  the  data-trend  of  Figure  2  but  is  too 
low  by  roughly  a  factor  of  two.  More  exact  analysis,  using  computer  codes 
for  calculating  propagation  loss  and  adding  the  reflection  loss  per  bounce, 
show  that  the  approximation  is  quite  accurate. 

The  fact  that  the  Marsh  formula  involves  the  correlation  length  implies 
that  the  scatter  beam-pattern  is  important  in  this  approximation  regime. 
For  kL<2n'  the  beamwidth  of  the  scattered  energy  tends  to  become  large. 
Much  of  the  energy  will  then  escape  from  the  channel  at  each  incidence  and 
be  absorbed  by  the  bottom. 

For  kL>2fl,  the  scattered  energy  approaches  specular  and  loss  goes  as  k2, 
which  is  characteristic  of  theEckart  [7]  regime.  While  perturbation  theory 
requires  kh  sin8  «1,  the  Eck8rt  solution  is  not  limited  by  this  criterion 
8nd  the  loss  approximation  becomes  "exact". 

In  the  low  frequency  limit  where  kL«2Tf,  scattering  becomes  diffuse  arid 
the  loss  goes  as  k4.  This  is  characteristic  of  Rayleigh  scattering  where  all 
scales  become  small  compared  to  acoustic  wavelength.  For  the  correlation 
length  44m,  rolloff  becomes  important  below  about  50  Hz. 

A  formula  based  on  the  method  of  small  perturbation  (MSP)  and  valid  in 
all  domains  will  be  used.  Equation  9.6.3  of  Brekhovskikh  and  Lysariov  [8] 
calculates  scattering  loss  by  a  2-dimensiorial  integral  over  the  surface 
spectrum.  In  Appendix  A,  we  derive  an  equivalent  integral  over  scattering 
angles,  which  shows  where  the  energy  goes.  Using  the  analytic  spectrum 
and  assuming  isotropicity  reduces  the  computational  problem  to  a  single 
integral  over  depression  angle.  Numerical  evaluation  of  the  B&L  integral 
has  been  found  to  give  identical  results 

The  coherent  reflection  coefficient  Q  is  defined  here  8S  the  fractional 
intensity  of  the  reflected  wave  in  the  specular  direction.  It  is  assumed, 
temporarily,  that  the  interface  acts  like  a  rough  free-surface.  Assuming 
small  displacements  compared  to  the  vertical  component  of  the  acoustic 
wavelength,  MSP  theory  yields  the  following  result: 
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4  TSf  2  2  tf  oo 

SL=  k  sin0  JdG’sin  0'cos0'  Jd0|rdr^(r)J_(Qr)/27r 

O  -1  o 

£2=k  Icos28+cos20'-2cosf  cosO  cosO'l1  /2 

k=2flf/C  V(r)=h2  Q0r  Ki(ft0r) 

OO  y 

J  rdr¥(r)Jo(Qr)=2h  q2/(q2  +  q2)2  =s2 

O 

First  Integral 

2  J*d0  S  2  /TJ  =  4  h2(l  +  u2)/^  (  ( I +U  )2-V2)3/2=  s; 

U=(k/Qo)2(cOS20+COS20)  V=2(k/£2o)2COS0COS0' 

Second  Integral 

4  */2  2 

SL=k  sin0[d0’sin  0'cos0'  s\ 

o 

Final  numerical  integral 

Eckart  SL=(2kh  sin0)2  Ru»1  kL»l 

2  1/2 

Marsh  SL=3.3(kh)  sin0/(kL)  Ru«l  kL»l 

Rayleigh  SL=(2/3)(kh)2(kL)2sin0  kL«l 
L=2/Q„  Ru=kL  sin2(0/2) 

Asymptotic  Approximations 

The  theory  is  based  on  Bragg-scatter  matching  between  acoustic  (k)  and 
surface  (ft)  wavenumbers  as  a  function  of  the  incident  grazing-angle  8,  the 
scatter  grazing-angle  9'  and  the  azimuthal  angle  0  relative  to  the  direction 
of  the  incident  plane-wave. 

The  first  (radial)  integral  yields  the  2-dimensional  analytic  spectrum. 
The  second  (azimuthal)  integral  yields  a  function  similar  to  the  analytic 
1 -dimensional  spectrum.  This  leaves  only  a  single  integral  over  scatter 
angle  8'  to  evaluate  numerically.  Asymptotic  solutions  provide  a  check  on 
results.  The  parameter  Ru  is  8  measure  of  forward-scattering  directivity 
and  occurs  in  Rutherford  scattering. 
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Figure  9:  Ice  spectrum  integrands  vs  scatter  angle. 

The  curves  of  Figure  9  show  the  scatter-loss  integrand's  vs  the  scatter 
grazing-angle  in  the  frequency  range  10-1000  Hz  for  an  incident  grazing 
angle  10°.  Results  have  been  normalized  by  dividing  by  (kh)2  sin8  in  order 
to  limit  the  scale.  Scattering  is  seen  to  become  diffuse  at  low  frequencies 
and  to  approach  specular  at  high  frequencies. 

The  integrands  can  be  considered  as  scatter  from  an  area  of  radius  equal 
to  the  correlation  length  and  “coherent"  with  the  incident  plane-wave.  The 
summation  over  a  random  ensemble  of  areas  makes  the  average  scattered 
energy  incoherent. 

The  integrands  resemble  beam  patterns;  however,  they  are  not  the  usual 
cross-sections  but  are  8zimuthally-integrated  values  and  multiplication 
by  cosS*  is  required  for  flux  integration.  At  low  frequency,  the  integrand 
approaches  the  vertical-dipole  pattern  sin28'  multiplied  by  cose'. 

From  the  integrands,  one  can  estimate  the  fraction  of  scattered  energy 
remaining  in  the  channel.  For  a  critical  trapping  angle  of  10°,  for  example, 
more  than  1/2  would  be  lost  at  each  reflection  at  frequencies  below  1  kHz. 
After  multiple  reflections,  incoherent  levels  should  then  converge  rapidly 
to  negligible  values,  providing  that  the  bottom  is  sufficiently  lossy. 
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Figure  10:  Ice  spectrum  scattering  coefficients  vs  frequency. 

In  Figure  10,  the  values  obtained  by  numerical  integration  are  plotted  vs 
frequency  for  selected  grazing  angles  between  4  and  16°,  which  represents 
the  domain  of  interest  for  long-range  propagation.  Values  have  again  been 
normalized  by  dividing  by  (kh)2sin8,  making  both  the  Marsh  8nd  Rayleigh 
asymptotes  independent  of  grazing  angle.  The  dashed  lines  at  the  right  8re 
the  Eckart  asymptotes. 

Note  that  the  Marsh  curve  tends  to  become  inaccurate  below  100  Hz  in 
the  Rayleigh  regime  8nd  also  at  high  frequency  in  the  Eck8rt  regime. 

Loss  estimates  below  100  Hz  have  already  proven  too  low  by  more  than  a 
factor  of  two  in  the  Marsh  approximation.  The  integral  values  become  even 
smaller  with  decreasing  frequency.  Since  the  conditions  for  MSP  theory 
are  evidently  not  exceeded  in  this  range,  there  is  little  reason  to  suspect 
errors  in  the  approximations. 

Failure  to  account  for  attenuation  indicates  that  impedance  effects  must 
be  involved.  Dissipation  within  the  ice  layer  is  a  potential  loss  mechanism. 
However,  much  higher  backscatter  levels  are  also  observed  experimentally 
th8n  predicted  by  free-surface  theory,  which  suggests  that  the  effects  of 
impedance  on  scattering  must  be  more  significant. 
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Figure  1 1:  Interface  impedance  model. 

The  impedance  problem  is  illustrated  in  Figure  1 1  where  g  and  p  are  the 
shear  and  compressional  grazing-angles  in  ice,  respectively.  Neglecting  air 
impedance  and  dissipation,  reflection  is  total.  Dissipation  is  included  by 
taking  cs  and  Cjto  be  complex,  as  described  in  Appendix  B. 

The  plane-wave  reflection  coefficient  then  becomes  [21: 

Rz=  [-1+  sin6  (Zr+iZi)/pc]/[  1+  sinB  (Zr+iZi)/pc] 

In  order  to  get  any  scattering  by  interfacial  roughness,  there  must  be  an 
impedance  discontinuity  at  the  boundary.  It  is  instructive  to  normalize 
relative  to  Zo.  Consider  the  lossless  case.  When  Zi/Zo-X),  there  is  perfect 
reflection  with  a  phase  reversal  and  the  interface  becomes  effectively 
pressure-release.  For  Zi/Zo^/p,  the  interface  discontinuity  is  inertial 
and  vanishes  for  equal  densities  in  the  two  media. 


Figure  12  Impedance  ratio  Zi/Zo  vs  grazing  angle 

Figure  12  shows  Zi/Zo  as  a  function  of  grazing  8ngle  for  the  frequency 
range  50  to  300  Hz  in  50  Hz  steps.  Compressional  and  shear  speeds  in  ice 
have  been  taken  as  3000  and  1800  m/s,  respectively,  with  the  shear  loss 
tani=.02.  For  simplicity,  the  same  value  is  used  for  compressional  loss, 
but  effects  are  minor. 

At  low  frequencies,  Zi/Zo->p1  /p  for  all  angles.  As  frequency  increases, 

resonance  effects  of  the  compressional  and  shear  waves  become  evident 
around  the  critical  angles  60°  and  33°,  respectively.  For  grazing  angles 
greater  than  60°,  compressional  effects  dominate  and  reactanre  oscillates 
between  inertial  (+)  and  stiffness  (-)  control.  Similar  effects  are  seen 
between  30°  and  60°  where  shear  effects  dominate.  Below  30°,  inertial 
effects  dominate,  diminishing  to  zero  as  the  shear  resonance  begins  to 
develop. 

In  the  Arctic  sound-channel  propagation  problem,  grazing  angles  less 
than  15°  and  frequencies  below  1  kHz  are  of  interest  at  longer  ranges.  In 
this  domain,  there  is  o  simple  transition  from  the  mass-looded  condition, 
with  negligible  boundary  discontinuity,  to  an  effectively  pressure-release 
condition  with  perfect  reflection  at  the  interface.  The  problem  is  that  the 
transition  occurs  at  such  high  frequencies. 
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Figure  13:  Impedance  ratio  Zr/Zo  vs  grazing  angle. 

Figure  13  shows  curves  of  the  impedance  ratio  Zr/Zo  for  tanA=.02.  It  is 
clear  that  internal  losses  are  small  except  at  resonances,  where  multiple 
internal  reflections  occur.  Mass  loading  become  dominant  at  small  grazing 
angles  and  the  elastic  properties  vanish.  In  this  regime,  dissipation  losses 
become  quite  small. 

The  ratio  Zi/Zo  is  most  important  in  the  case  of  vertical-dipole  scatter 
because  it  is  a  measure  of  the  phase  shift  of  the  reflected  wave  relative 
to  the  ice/water  interface.  No  phase  shift  means  perfect  reflection  at  the 
interface  while  a  delay  equal  to  that  in  a  water  layer  of  equal  thickness- 
means  that  the  ice/water  impedance  discontinuity  vanishes.  Corresponding 
effects  in  the  time  domain  are  shown  in  Appendix  C. 

Perturbation  analysis  for  finite  impedance  19]  is  outlined  in  Appendix  D. 
Results  show  that  the  vertical-dipole  component  of  scatter  loss  Is  simply 
the  free-surface  integral  with  the  integrand  multiplied  by  the  effective 
reflection  coefficient  Ro2. 

The  analysis  also  shows  that  there  is  second  vertical-dipole  component 
of  the  velocity  type,  analogous  to  a  piston  in  a  baffle.  The  scatter  loss  for 
the  infinite-impedance  case  involves  a  similar  integral.  For  the  finite- 
impedance  boundary,  the  integrand  is  simply  multiplied  by  the  effective 
reflection  coefficient  R»2 
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Verticol-dipole  models 


Figure  14:  Attenu8tion-vertical  dipole  models. 

Figure  14  compares  predictions  of  the  various  vertical-dipole  models  to 
the  attenuation  data.  The  gradient  is  taken  as  0.06/s  and  the  grazing  angle 
as  8°,  making  the  skip-distance  roughly  7  km. 

The  dashed  line  is  the  Marsh  formula.  Note  that  the  predictions  are  low 
by  more  than  a  factor  of  two  at  low  frequencies. 

Curve  I  is  the  numerical-integration  result  for  the  free-surface  model. 
The  slow  low-frequency  rolloff  involves  transition  to  the  Rayleigh  regime 
where  scattering  becomes  diffuse.  A  slow  transition  to  the  Eckart  regime 
of  "specular  scatter  is  seen  at  the  higher  frequencies. 

Curve  V  shows  the  effects  of  the  impedance  factor  Ro2  on  curve  1.  Note 
that  the  attenuation  decreases  rapidly  in  the  mass-loading  regime  below 
about  200  H2.  The  “piston"  component  turns  out  to  be  negligible  throughout 
the  entire  range  and  is  omitted. 

Dissipation  loss  is  calculated  from  the  magnitude  iRzl2  Curve  A  shows 
the  results  for  tanA=.02.  Total  attenuation  is  then  the  sum  of  A  and  V. 

Clearly,  the  results  do  not  agree  with  the  data  and  “improvements”  seem 
to  have  made  matters  much  worse.  Analysis  of  seismometer  data  from 
experiments  by  Buck  [101  suggests  8  possible  reason  for  the  failure. 
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Seismometer  experiments 


Perfect  acoustic  contact  8t  the  ice/water  interface  has  been  assumed. 
Trapped  air  could  make  the  impedance  very  low.  Comparison  of  hydrophone 
and  seismometer  signals  shows  that  this  is  not  the  case. 

Figure  15  shows  the  response  (greatly  exaggerated)  of  the  ice  layer  to  an 
incident  pl8ne-wave.  The  transfer  function  of  pressure  to  surface  velocity 
can  be  estimated  by  comparing  the  signals  received  from  a  distant  source. 
At  low  frequencies  and  grazing  angles,  thin-plate  theory  is  valid  and  the 
motion  of  the  median  plane  (dashed  line)  is  vertical.  Neglecting  the  small 
difference  in  densities,  the  vertical  component  of  surface  velocity  is  the 
same  as  the  free-surface  case:  Uv~(2P/pc)  sin8,  where  P  is  the  pressure. 
The  horizontal  component  is  Uhs  iUv(kD/2)cos8,  where  D  is  the  thickness. 

Figure  16  compares  hydrophone  and  vertical  seismometer  signal  levels 
at  25  and  50  Hz.  The  equivalent  plane-wave  pressure  is  estimated  from  the 
hydrophone  response  at  the  experimental  depth  100m.  Agreement  with  the 
estimated  response  for  8=8°  (dashed  line)  is  seen  to  be  reasonable. 

Figure  17  compares  the  horizontal  and  vertical-seismometer  signals. 
Estimates  for  the  experimental  condition  D=5m  are  shown  by  the  dashed 
lines.  The  fact  that  horizontal  velocities  are  much  greater  than  predicted 
indicates  a  “rocking"  motion  involving  slopes.  The  associated  horizontal- 
dipole  term  of  Appendix  D  must  therefore  be  considered. 
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Figure  16:  Vertical  seismometer  vs  hydrophone  data. 
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Figure  17:  Horizontal  vs  vertical  seismometer  data. 
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Figure  16:  Attenuation-general  impedance  model. 

The  scatter-loss  integral  for  the  horizontal-dipole  case  in  Appendix  D  is: 

SU2  k2  cos8  cotB  0(€/2d8*  cos8’  Roo2  _  J*d$  Qx2  S2(fl)/rt 

where  Roo2  is  the  effective  reflection  coefficient  for  the  rigid-boundary 
solution  and  Ox=k  (cos8-cos8'cos$)  is  the  Bragg-scattering  condition 
associated  with  the  x  component  of  the  slopes.  Since  an  analytic  solution 
of  the  $  integral  has  not  been  found,  the  double  integral  must  be  calculated 
numerically.  Note  that  SL-»co  as  8-»0  only  for  the  rigid-bound8ry  case  and 
remains  finite  for  finite  impedance  since  Roo  contains  a  sin^  term. 

The  attenuation  curves  shown  in  Figure  18  are  calculated  for  the  same 
conditions  as  Figure  14.  Curve  H  is  the  horizontal-dipole  component,  curve 
V  is  the  vertical-dipole  component  and  curve  A  is  dissipation. 

It  is  clear  that  the  hori2ont8l-dipole  scatter  component  is  8  significant 
attenuation  mechanism  at  the  lower  frequencies.  However,  predictions  do 
not  agree  with  the  data  in  this  range.  Sensitivity  to  ray  angle  should  be 
noted.  Since  angle  effects  are  quite  complicated,  a  detailed  analysis  by 
either  norm8l-mode  or  numerical  methods  may  be  indicated.  Sensitivity  to 
the  PDF  of  draft  is  more  readily  examined. 
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Figure  19:  Composite  general-impedance  model. 

Figure  19  shows  that  predictions  can  be  put  into  only  slightly  better 
agreement  with  data  using  the  composite  PDF  model  of  Figure  8;  i.e.  a 
bimodal  model  with  standard  deviations  1.1m  and  3.3m,  weighting  factors 
0.7  and  0.3  and  mean  drafts  4m  and  8m,  respectively.  Values  for  the  two 
modes  were  calculated  numerically  and  the  area-weighted  losses  added. 

The  m8in  effect  of  the  composite  model  is  due  to  the  impedance  of  the 
dominant  8m  component  of  draft.  The  higher  impedance  causes  an  increase 
in  low-frequency  loss  and  a  reduction  at  high  frequencies,  resulting  in  a 
slightly  better  fit  to  the  data. 

Obviously,  the  agreement  can  be  improved  by  further  modification  of  the 
PDF  model;  e.g.  a  multimodal  model  can  be  devised  that  might  provide  a 
better  overall  fit.  However,  in  view  of  the  problems  involving  r8y-angle 
sensitivity,  as  well  as  questions  about  data  accuracy,  such  refinements  do 
not  appear  justified  at  this  time  nor  are  the  prospects  for  success  very 
encouraging. 
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Conclusion 


It  is  apparent  from  this  analysis  that,  although  scattering  at  the  rough 
ice/water  interface  remains  the  most  likely  mechanism  for  the  excess 
attenuation  noted  in  Arctic  propagation  experiments,  consideration  of 
finite-impedance  effects  does  not  appear  to  account  for  the  data. 

The  ice-reflectivity  problem  has  turned  out  to  be  quite  complex.  At 
finite-impedance  boundaries,  scattering  involves  horizontal-dipole  and 
'piston"  components  in  addition  to  the  usual  vertical-dipole  component 
associated  with  free  surfaces.  While  the  horizontal-dipole  may  be  the 
most  significant  loss  mechanism  at  lower  frequencies,  predicted  values 
are  still  much  too  low. 

The  scattering  theory  is  based  on  the  usual  Bragg  condition  wherein  the 
scatter  angles  are  determined  by  matching  of  the  acoustic  and  roughness 
wavenumbers.  Rigorously,  the  MSP  approximation  requires  small  values  of 
the  Rayleigh-roughness  parameter.  This  is  violated  mainly  in  the  Eckart 
regime,  where  the  solution  is  not  limited  by  the  criterion.  Approximation 
errors  are  not  expected  to  be  serious  much  below  1  kHz. 

Effects  of  finite  impedance  on  vertical-dipole  scattering  are  quite  clear 
at  lower  frequencies  and  grazing  angles  where  inertial  forces  dominate.  In 
this  regime,  scattering  arises  from  spatial  variations  in  mass.  Since  ice 
displaces  an  equal  (almost)  mass  of  water,  scattering  tends  to  vanish.  The 
scatter  loss  tdB/bounce)  is  approximately  proportional  to  grazing  angle. 
The  skip  distance  has  similar  dependence,  making  loss  in  dB/km  nearly 
independent  of  angle. 

In  the  horizontal-dipole  case,  the  dB/bounce  loss  is  a  more  complex 
function  of  grazing  angle,  making  the  analysis  of  propagation  much  more 
sensitive  to  angle  effects.  However,  since  predictions  are  far  too  small  in 
the  inertial  regime  below  about  100  Hz  for  all  reasonable  angles,  current 
prospects  for  an  adequate  predictive  model  appear  rather  remote. 
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Appendix  A:  Free-surface  theory 


Perturbation  theory  for  scattering  at  a  free  surface  is  treated  in  Ref.  8. 
The  surface  boundary  condition  is  given  by  Eq.  9.2.4  as: 

p+  £  dp/dz=0  1 A 

where  p  is  the  acoustic  pressure,  £  is  the  surface  displacement  and  z  is 
depth.  Let  p=  po+  ps  where  po  is  the  pressure  for  £=0  and  ps  is  the  scattered 

pressure  for£*0.  Equating  the  zero  and  first-order  perturbation  terms  as 
in  Eqs.  9.2.6-7,  we  have: 

P0=°  ps+  K  ap0/dZzO  2A 

The  scattered  field  corresponds  to  vertical-dipole  reradiation.  Instead 
of  following  the  derivation  of  Ref.  8,  we  take  the  more  physically  oriented 
approach  outlined  below. 

The  farfield  component  of  the  Green's  function  for  a  scattering  element 
and  surface  image  yields  the  differential  scattered  pressure: 

dPs*ikP0sine  dA  Iexp(ik<;  sin8')-exp(-ikf;  sine')]  exp(ikR)/2flR  3A 

where  k  is  the  acoustic  wavenumber,  Po  is  the  incident  pressure,  8  and  8' 

are  the  incident  and  scatter  grazing-angles,  dA  is  an  element  of  surface 
area  and  R  is  the  radius  to  a  distant  field  point.  For  small  values  of  the 
exponential  arguments,  Eq.  3A  becomes: 

dPss-Po  k2 £  sin8  sin8'  dA  exp(ikR)/TJR  4A 

For  a  plane  wave  incident  in  the  x  direction,  we  then  have: 

Ps~-P0  k2  JdxJdy  £(x,y)  sin8  sin8‘  exp[ik(R+x  cos8)]/rfR  5A 

where  x,y  are  the  surface  coordinates.  The  scattered  field  is  seen  to  be 
equivalent  to  broadside  (y  integral)  and  vertically-steered  (x  integral) 
reradiation  by  a  phased-array  of  vertical  dipoles. 

The  Bragg-scatter  condition  for  matching  the  wavenumber  spectrum  of 
£(x,y)  is  Q  =k(cos8-cos8‘  cos$)  and  Q=k  cos8‘  siruji.  Continuous  spectra 

x  y 

involve  mean-squared  values  and  spatial  convolution  of  Eq.  5A  is  required. 
For  the  case  of  isotropic  roughness,  it  can  be  shown  that  the  result  is: 

Ps2(8'i)~  P02  Ao  k4  sir^B  sinks’ 0(°°rdr  y(r)  J0(Or)  /(nR)2  6A 

where  q=(ox2+  oy2)1/2,  r=(x2+y2)1 12 ,  ¥(r)  is  the  correlation  function  of  £ , 

sin2  8‘  is  the  scattering  beam-pattern  of  an  element  dA,  Ao  is  the  effective 

"coherent"  area  and  Q=k  [cos28+cos28'-2  cos8‘  cos8  cosifrl1^2  is  the  Bragg 
wavenumber  condition. 


0*7 

/LO 


7A 


The  total  flux  scattered  by  the  area  is: 

Ps2Ao  R2  0J*/2  d8‘  C0S8'  _J*d$  Ps2(8'i) 

We  defined  the  intensity  reflection-coefficient  as  Q=1-SL,  where  SL  is 
scattering  loss.  Conservation  of  total  vertical  flux  8t  the  surface  requires 
P02(1-Q)  sin0=  Ps2.  Averaging  over  a  random  ensemble  of  areas  and  taking 

SL=<Ps2>/Po2  sin8  gives: 

SU2  k4  sinO  0|*/2  d0'  sin2©’  cosO'^j’d*  oJ°°rdr¥(r)  J0(Qr)/id  8A 

The  radial  integral  gives  the  surface  spectrum  S2(q)  and  Eq.  8A  becomes: 

SL~2  k4  sin8  0("/2  d8’  sirr^'  cos0'_J*dt  S2(ft)/fl  9A 

Equation  9A  is  equivalent  to  Eq.  9.6.3  of  Ref.  8.  In  terms  of  the  present 
model  parameters,  the  E&L  integral  is: 

SL~2k2  sin0  [q  dQ  )<$!>  [sin28+2(a/k)cos8  cos^-f^/k)2]1  /2  S2(Q)/rt  1 0A 

where  the  radical  derives  from  the  Bragg  condition  and  the  integration 
range  is  over  real  values.  The  first  two  terms  of  the  radical  correspond 
uniquely  to  the  Eckart  and  Marsh  regimes  while  the  last  term  causes  the 
low-frequency  rolloff  in  the  Rayleigh  regime. 

Equivalence  of  the  two  formulae  is  seen  by  changing  the  integration  of 
Eq.  10A  from  wavenumber  to  angle  variables  8'  and$  as  indicated  in  Ref.  8. 
From  Eq.  9.3.10,  Q  dQ  becomes  k2  sin8'  cos8'  d8'  and,  from  Eqs.  9.6. 1-3, 
[sin28+2(Q/k)cos8  cos^-fo/k)2]1  /2=sin8'. 
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Appendix  B:  Complex  impedance 


In  Figure  1 1,  tangential  forces  are  assumed  zero.  For  no  losses,  matching 
of  horizontal  wavenumbers  at  the  ice/water  interface  gives  the  relations: 
k  cos8  =  k^osp  =  k^cosg 

on  cosp  =  (c^c)  cos8  cosg  =  (cs/c)  cos8  IB 

Then:  sinp  =  11-  (ct  /c)2  cos28]1  /2 

sing  =[  1  -  (cs/c)2  cos20]1  i2  2B 

The  waves  in  ice  become  evanescent  for  imaginary  values  of  the  radicals. 

In  evaluating  the  equation  of  Figure  1 1,  we  make  use  of  the  identities: 
(sin2q)2  =  (2  sing  cosg)2  =  4  cos2g  ( 1  -  C0S2g) 

(cos2g)2  =  (2  cosg  - 1)2  3B 

The  impedance  equation  is  derived  from  Eq.  10.7  of  Ref.  2  by  taking  the 
impedance  of  the  air/ice  interface  to  be  zero.  Note  that  the  equation  has 
been  changed  to  N2  -  M2  in  place  of  M2  -  N2  in  order  to  correct  for  a  change 
in  sign  of  the  arguments.  That  this  is  correct  is  seen  from  the  asymptotic 
values  Z=i  p1c1tan(k1D)  at  normal  incidence  and  ZH>i  p^D,  cdh>0,  where  the 

latter  is  the  mass-loading  limit  8nd  pjD  is  mass  per  unit  area. 

Losses  are  included  by  substituting  c1  ( 1  +i )  and  cs(1+iAs).  Values  of  a 

are  both  taken  to  be  .02  for  simplicity.  The  BASIC  program  for  calculating 
Zi/Zo  and  Zr/Zo  is  shown  in  Figure  IB  The  functions  used  in  the  subroutine 
are  sinx/x  and  tanx/x,  where  x  is  the  appropriate  shear  or  cornpressional 
argument.  The  corresponding  coefficients  in  the  main  program  have  been 
changed  accordingly. 

In  all  calculations,  impedance  is  approximated  as  that  of  the  unperturbed 
boundary. The  method  of  Lapin  [IB]  is  more  general;  however,  resonance 
angles  occur  where  the  wavenumbers  of  the  scattered  field  match  the 
shear  and  cornpressional  modes  in  the  ice  layer.  Numerical  solutions  of  the 
latter  equations  indicate  that  the  errors  are  not  serious,  particularly  at 
lower  frequencies  where  impedance  is  nearly  independent  of  angle. 

Ref.  IB:  A.  D.  Lapin,  "Scattering  of  sound  by  a  solid  layer  with  rough 

boundaries",  Sov.  Phys.  Acoust  J_2  46-51  (1966) 
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'Ice  impedance- Brekhovskikh 

2000  'Subroutine 

dG=4‘  ice  thickness 

kkr=k/r/(1+dA2) 

03=  1 800  ’shear  speed 

kki=-kkr*d 

cc=3OO0  ’compressional  speed 

pp=1 -(r*c)A2*( 1-dA2) 

ds=.02 loss  tan  shear 

qq=-2*d*(r*c)"2 

dc=ds  loss  tan  compressional 

v=(ppA2+qqA2)A(  1  /4) 

rs=C3/1 500:rc=cc/1 500 

ii=ATN(qq/pp) 

xscale=3:yscale=  1 

IF  pp>0  THEN  2100 

GOSUB  1000 ‘Graphics 

ii=ii+3. 141592654* 

theta=8 

2100 

t=theta/57.3:s=SIN(t):c=C0S(t) 

sr=v*C0S(ii/2> 

aO=2*(rs*c)A2 

si  =w*SIN(ii/2) 

br=2*rs*a0*(  1  -  3*dsA2)  'shear  coefficients 

kr=kkr*sr-kki*si 

bi=2*rs*aO*ds*(3-dsA2) 

ki=kki*sr  +  kkr*si 

al =30*( 1 -dsA2)- 1 

cu=C0S(kr) 

ur=a1A2-  4*aOA2*dsA2 

su=SIN(kr) 

ui=4*a0*a1*da 

jj=EXP(ki) 

a2=(rc*c)A2 

shu=(jj-1/jj)/2 

u=1-a2*(1-dcA2) 

chu=(jj+  1/jj)/2 

v=-2*a2*dc 

w=(su*chu)A2+(cu*shu)A2 

v=u*2+vA2 

ur=3U*cu/w 

vr=(u+v*dc)/w 

ui=-shu*chu/v 

vi=(-v+u*dc)/v 

vr=3u*chu/v 

ar=rc*(ur*vr-ui*vi)  'compressional  coefficients 

vi=-cu*3hu/v 

ai=rc*(ui*Yr+vi*ur) 

Mr=ur*$r-ui*si 

FOR  n=OTO  3  STEP  .1 

Mi=ur*si  +  ui*sr 

f=50*(n+1)  ‘acoustic  frequency 

Nr=vr*sr-vi*si 

k=2*3.1 4*f*d0/1 500  Vavenumber*dO 

Ni=vr*si+vi*sr 

r=rs:d=ds 

GOSUB  2000 

Msr=Mr.Msi=Mi:Nsr=Nr:N3i=Ni 

r=rc:d=dc 

GOSUB  2000 

Mcr=Mr:Mci=Mi:Ncr=Nr:Nci=Ni 

Mr=ar*Mcr-ai*Mci  +  br*Msr-bi*M3i 
Mi=ar*Mci+ai*Mcr+br*M$i  +  bi*Msr 

Nr=ar*Ncr-8i*Nci  +  br*Nsr-bi*Nsi 
Ni=ar*Nci+ai*Ncr+br*Nsi+bi*Nsr 
zi=Mr-(Mr*(NrA2-NiA2)  +  2*Mi*Nr*Ni)/(MrA2+MiA2) 
2r=Mi  +  (Mi*(NrA2-NiA2)-2*Mr*Nr*Ni)/(f1rA2+MiA2) 
2i»-2i*.9/k:2r=2r*.9/k  ‘normali2e//20 

x=xO+n*1 

yi=yO-zi*h:yr=yQ-zr*h 

LINE  (xt,ytr)-(x,yr):LINE  (xt,yti)-(x,yi) 
xt=x:utr=yr:yti=yi 

NEXT  n 

STOP 

RETURN 

Figure  IB:  BASIC  program  for  Zi/Zo  and  Zr/Zo. 
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Appendix  C:  Impulse  response 


Figure  1C  shows  the  impulse  response  of  the  ice  layer  as  a  function  of 
grazing  angle.  The  waveforms  8re  calculated  by  inverse  FFT  (N=256)  of  the 
reflection  coefficient  Rz. 

A  gaussiari  waveform  is  used  for  the  incident  impulse  (not  shown),  which 
is  done  by  gaussian  weighting  the  Fourier  amplitudes.  This  serves  to  limit 
the  peak  amplitude  and  an  also  to  prevent  ringing  at  high  frequency.  The 
non-dimensional  Fourier  frequency  of  the  fundamental  is  taken  as  FD=1. 
Overlap  effects  of  adjacent  periods  appear  negligible. 

A  time  shift  is  needed  in  order  to  show  the  precursor  properly.  Instead 
of  changing  the  impulse  time,  we  take  n=50  as  the  start  time  and  make  use 
of  FFT  periodicity. 

The  parameters  used  are  the  same  as  those  in  Figs.  12  and  13.  The  shear 
critical  angle  is  roughly  33°.  Individual  reflections  at  both  boundaries  are 
seen  at  the  larger  grazing  angles.  The  impulses  merge  at  the  critical  angle 
where  the  shear  wave  becomes  evanescent.  Polarity  of  the  high-frequency 
components  then  becomes  negative.  Inertial  effects  at  lower  frequencies 
and  grazing  angles  is  evidenced  by  the  persisting  slower  oscillations. 
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Appendix  D:  Finite-impedance  theory 


n . i“4  1 1  l—  i— . 1 i-v  +  u  ^ .  I- 1  ■  +" k .  +  L-.  f  «.— j  L.— i  ,-i  * i  .  t  H-+  ft  -i +  I- ft  +  ,-i  » »ft  n ♦  o 

rei  uji  uauuii  uieui  y  i  ui  uuuii  i  ee  aiiu  i  ty  iu  sui  i  aues  is  u  eaieu  m  r.ei .  u. 


"he  respective  boundary  conditions  are: 


p+  <  dp/dz=0 

dp/dz+  £  d2p/dz2  -  7p  ■  %  -  0  ID 

where  7=i d/bx  +  jd/dy.  In  the  rigid-boundary  case,  the  slopes  %  involve 
horizontal -dipole  reradiation  while  the  displacements  £  involve  vertical 
velocity-dipole  reradiation  analogous  to  a  piston  in  a  baffle.  Equating  the 
terms  of  the  two  perturbation  orders  to  zero  gives: 

Po=0  Ps+  y  3po/dZ=G 


dpJdZ-0  Gps/3z+  £  d2po/dZ2  -  Vpo  ■  7y  =  0  2D 

Theory  for  t'inite-imipedance  boundaries  is  treated  in  Ref.  9.  The  analysis 
involves  normal-modes  and  is  of  limited  use  in  the  plane-wave  reflection 
case;  however,  the  finite-impedance  perturbation  equation  is  given  as: 

p+  £  dp/dz=  (dp/dz+  <;  d2p/dz2  -  7p  •  7<)  (Z+  f  dZ/dD)  3D 

where  Z~Zi/wp  (Zr  can  be  neglected).  Equations  analogous  to  Eq.  2D  become: 
po-  zdpo/dz=0 


pf-  Zdp./dz  =  y  [(- 1  +  dZ/dD )dp./r)z+  Zd2p./dz2]-  Z7po  •  7 y  4D 

The  vertical-dipole  component  becomes. 

pc(1-ikZ  sin0')«  y  (-1+dZ/dD)  dpjdz-  -  y  [1-dZi/dZo]  dpo/dz  5D 

which  is  the  free-surface  result  with  impedance  correction  factors.  From 
the  relations  kZ=Zi/pc,  po=Po(1+Rz),  dpo/dz=ikP0(1-Rz)  sin9,  the  intensity 

coefficient  for  finite  Z  becomes: 


Ro2=[1-dZi/dZoF/{(1+(sin0'  Zi'/pc)2]  M+(sin8  Zi/pc)2]}  6D 

which  becomes  a  reflection  factor  in  the  free-surface  integrand.  Errors  at 
lower  frequencies  are  evidently  not  serious  if  we  take  Zi'=Zi. 

The  velocity-dipole  term  Z(dps/dz+yd2po/dz2)  is  orthogonal  and  can  be 


treated  independently.  By  the  same  method,  the  solution  for  Z=oo  is: 
Pc2(0\'t>)~  P  2  k4  A  si n4G  S0(Q)/(dR)2 


7D 


where  the  element  beam-pattern  is  now  semi -isotropic.  Averaged  over  a 
random  ensemble  of  areas,  the  scatter-loss  integral  for  finite  z  becomes: 
SU2  k4  sin38  .  |ff/2d0'  cos8'  R<x>2  bd<|>  S0(Q)/rf 

V  f  "IT-  A. 


8D 


=k4  sin3 8  0\wUm‘  cos0'  Roo?  si  '(Q) 

where  Roo2=(sin8'  Zi'/pc)2  (sine  Zi/pc)2/{[  1 +(sin0‘  Zi7pc)2j[  1  +(sin8  Zi/pc)2]} 
is  the  effective  reflection  factor. 
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By  the  same  method,  the  horizontal-dipole  solution  forz=«»  is: 

Ps2(8',<j>)~  Po2  k2  A0  cos20  Qx2S2(Q)/(tTR)2  9D 

where  Qx=k  (cos8-cos8'cos4i)  is  the  Bragg  condition  associated  with  the 
gradient  of 

Results  are  quite  different  from  the  vertical-dipole  case.  At  the  lower 
frequencies,  most  of  the  scatter  is  backwards.  Forward  scatter  increases 
with  increasing  frequency  and  backscatter  decreases.  The  element  beam 
pattern  is  not  the  usual  horizontal-dipole  pattern  because  the  excitation 
phase  is  delayed  along  the  x  axis  and  the  dipole  null  occurs  in  a  semi-cone 
(ftx=0)  including  the  specular  direction. 

Averaged  as  in  Appendix  A,  the  scatter-loss  integral  for  the  horizontal- 
dipole  component  with  finite  Z  is: 

SLs2  k2  cos8  cot8  oJ^2d0’  cos8'  Roo2  _  J"d4i  f2x2  S2(q)/tT  1 0D 

Assuming  no  coherence  between  slopes  and  displacements,  this  component 
is  independent  and  additive  also. 

Equivalent  equations  for  the  velocity  8nd  horizontal-dipole  components 
can  also  be  derived  by  the  method  of  6&L  [8].  Forz=oo  the  equation  for  the 
velocity-dipole  component  is: 

SLs2  k2  sin38  Jq  dQ  Jd$  S2fo)/J7sin8'  1 1 D 

and  the  equation  for  the  horizontal-dipole  component  is: 

SL~2  cos8  cot8  }  Q3dft  {tty  cos2!})  S2(fi)/Tfsin8‘  1 2B 

where  sin0’=[sin28+2(o/k)cos8  cos^-fo/k)2]1 /2  and  the  range  of  integration 
is  over  real  values  of  the  radical. 

The  rigid-boundary  case  is  also  treated  in  Ref. ID  by  a  different  method. 
The  asymptote  in  the  Eckart  regime  is  shown  to  be  SL=(2kh  sin8)2,  which 
is  the  same  (magnitude)  as  the  free-surface  case.  This  result  is  obtained 
from  Eq.  1 1D,  keeping  only  the  first  term  in  the  radical.  The  asymptote  in 
the  Marsh  regime  is  SL « (kh)2/(kL)3/2  sin8.  This  result  is  obtained  from 
Eq.  1 2D,  keeping  only  the  second  term  in  the  radical.  The  asymptote  in  the 
Rayleigh  regime  is  SL  c*  cos8  cot8  (kh^/tkL).  This  result  is  also  obtained 
from  Eq.  12D  by  including  the  third  term  in  the  radical.  Numerical  results 
of  Eqs.  8D  and  10D  for  Zi=oo  have  been  found  to  be  identical. 

Ref. ID:  V.  P.  Lysanov,  “Scattering  of  Sound  by  Irregular  Surfaces",  in 

Ocean  Acoustics  ed.  L.  M.  Brekhovskikh  (Moscow,  1974)  Part  IV. 
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